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Abstract 

. Special features of magnetohydrodynamic waves linear dynamics in smooth shear flows are stud- 

P 

led. Quantitative asymptotic and numerical analysis are performed for wide range of system 
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> 

X 

■ magnetosonic and Alfven waves at high shear rates. Studied phenomena should be decisive in the 

elaboration of the self-sustaining model of magnetohydrodynamic turbulence in the shear flows. 
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parameters when basic flow has constant shear of velocity and uniform magnetic field is parallel 
to the basic flow. The special features consist of magnetohydrodynamic wave mutual transfor- 
mation and over-reflection phenomena. The transformation takes place for arbitrary shear rates 
and involves all magnetohydrodynamic wave modes. While the over-reflection occurs only for slow 
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I. INTRODUCTION 



Flows with inhomogeneous velocity profiles are one of the prevalent and still not com- 
pletely understood macro-systems. They occur in laboratory experiments, industrial appli- 
cations, the earth's atmosphere, oceans and many astrophysical objects. We concentrate 
attention on smooth shear flows (flows without inflection point in the velocity profile) that 
are linearly stable according to the canonical hydrodynamics. Faults in the comprehension 
of shear flow dynamics is caused by the incompleteness of the classical method of normal 
modes. Specifically, method of modes leaves out of account the non-normality of the linear 
operators that often produces powerful transient development of a subset of perturbations 
P, 0, 0, y, 0, ^. Indeed, in the case of the non- normality of inear operators, correspond- 
ing eigenf unctions are not orthogonal and strongly interfere js]. Hence, a correct approach 
should fully analyze the eigenfunction interference, which is feasible in asymptotic: in fact, 
in modal analysis only the asymptotic stability of flow is studied, while no attention is paid 
to any particular initial value or finite time period of the dynamics - this period of the 
evolution is thought to have no significance and is left to speculation. 

Recognition of the importance of the non-normality in the linear stability resulted im- 
pressive progress in the understanding of the shear flow phenomena in 90s of the last century. 
The early transient period for the perturbations has been shown to reveal rich and compli- 
cated behavior leading to various consequences. It has been grasped phenomena that are 
overlooked in the framework of the modal analysis. 

The progress in the understanding of the shear flow phenomena has been achieved us- 
,„g so-caUed „o.-,„oda. approach (see. e.g.. Refs. flft). wMch implicates a change of 
independent variables from a laboratory to a moving frame and the study of temporal evo- 
lution of spatial Fourier harmonics (SFHs) of perturbations without any spectral expansion 
in time. Strictly speaking the non-modal approach is applicable to the flows with constant 
velocity shear. However, it is obvious that the results obtained in the framework of the 
approach are valid for any shear flows without inflection point, if the wavelength of modes 
is appreciably shorter than the characteristic length scale of the inhomogeneity. 

The progress involves novel results on the time evolution of vortex and wave perturba- 
tions; elaboration of a concept of self-sustainingturbulence in the spectrally stable shear 
flows (that was labelled as bypass concept) [tI Isl lol llf| . Among the novelties one should 
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stress (in the context of the present paper) disclosure of the hnear mechanisms of reso- 

nnn 

nance transformation of waves at low shear rates |^ I12I. 13| and non-resonant conversion of 
vortices to waves at moderate and hieh shear rates [14 L The later phenomena should be 
inhere, .o How .y.ems. Th. was shown in pape.^ Q. wKe.e appearance, of so,ne 
waves in solar wind and galactic disks are explained by the linear wave transformations in 
the magnetohydrodynamic (MHD) flows. 

Aim of the present work is quantitative and qualitative study of the linear evolution of 
three dimensional MHD waves at low and high shear rates. Equilibrium density, pressure 
and magnetic field are assumed to be homogeneous. Specifically, we perform analytical and 
numerical study of the important effect caused by linear forces - mutual transformations of 
MHD waves. The waves are coupled (during a limited time interval) at the specific system 
parameters and mutual transformations of two, or even all three of them occur. Analytical 
expressions of transformation coefficients are obtained for wide range of the system parame- 
ters. We study the wave over-refiection phenomena at high shear rates: MHD waves extract 
flow energy, are amplifled non-exponentially and are over-reflected. 

Mathematical methods used in this paper are similar to ones that were flrst developed 
for quantum mechanical problems: non-elastic atomic collisions [1™ and non-adiabatic tran- 

nn 

sitions in two level quantum systems jl7|, 118||. Later, the same asymptotic methods were 
successfully apphed to various problems (see, e.g., Ref. 1^) including interaction of plasma 
waves in the media with inhomogeneous magnetic field and/or density j^]. 

The paper is organized as follows: employed mathematical formalism is presented in Sec. 
II. General properties of resonant transformation of wave modes that takes place at small 
shear rates, are presented in Sec. HI. Particular cases of the resonant transformation of 
MHD wave modes are described in Sec. IV. In Sec. V dynamics of SFHs of the wave modes 
at high shear rates and over-refiection phenomenon is studied. Conclusions are given in Sec. 
VI. 



II. MATHEMATICAL FORMALISM 



Consider compressible unbounded shear fiow with constant shear parameter 
(Uo(A?/, 0, 0)) and uniform density (po), pressure (Pq) and magnetic field directed along 
the streamlines (Bo||Uo). Linearized ideal MHD equations governing the evolution of den- 
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sity (p'), pressure {p'), velocity (u') and magnetic field (b') perturbations in the flow are: 



(5t + Uo- V)p' + poV-u' = 0, (1) 



{dt + Uo ■ V)u' + (u' • V)Uo + -Vp' = -^Bo X (V X b'), (2) 

Po 47rpo 



(dt + Uo • V)b' - (Bo • V)u' + Bo(V ■ u') = 0, (3) 

V ■ b' = 0. (4) 

Further we use the standard technique of the non- modal approach i.e., expand the 
perturbed quantities as: 

t) = ip'{k, t) exp [i{Kx + ky{t)y + k^z)] , (5) 

where ijj' = {u^' ,Uy\uJ , p\p\hj ,hy\hj), ky{t) = ky{0) - k^At, and {kn,, ky{0), k^) are 
the wave numbers of SFHs at the initial moment of time. It follows from Eq. (jSJ, that 
ky{t) varies in time. This fact can be interpreted as a "drift" of SFH in phase space. This 
circumstance is caused by the fact, that perturbations cannot have a simple plane wave form 
in the shear flow due to the shearing background j^. 

Using the following thermodynamic relation p' = Cg'^p', where Cg is sound speed, and 
introducing non-dimensional parameters and variables: 

S = A/{VAk,), T = kM, (i = cl/Vl 
Kz = kz/kx, Kyi^T) = ky/kx St, 

(6) 

K\r) = 1 + K'yir) + Kl p(k, r) ^ ^(k, r)/po, 
b(k, r) = zb'(k, r)/Bo, v(k, r) = u'(k, t)/Va 

(where Va = -Bo/V^^rpo is the Alfven velocity) the set of Eqs. ©-(jll) is reduced to: 

^ = + Ky{T)vy + KzVz, (7) 



^ = -Svy - (3p, (8) 
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dr 



-(3Ky{T)p + (1 + K'^{T))by + Ky{T)KA, (9) 



cil) 

= -(3K,P + (1 + Kl)h, + Ky{T)K,hy, (10) 



d7 = (11) 



K + Ky{T)hy+K,h,=Q. (13) 

SFH energy density in the non-dimensional form may be defined as a sum of non- 
dimensional kinetic, magnetic and compression energy densities: 

E(k, r) ^ E\\i, r) + E'"(k, r) + E^(k, r), (14) 

where: 

3 

i?'=(k,r)^^K|V2, (15) 

i=l 

E-(k,r) = X^|6,|V2, (16) 
1=1 

i?'^(k,r)=/3|p|V2, (17) 

For further analysis it is convenient to introduce a new variable d = p + Ky{T)hy + KJ)z 
and rewrite the set of Eqs. (f7|)-()13|) in the following form [lj||: 



dr2 



+ Cnrf + Ci2(r)6, + Ci36. = 0, (18) 



dr2 



+ C^22(r)6, + C2i(r)d + C23(r)6, = 0, (19) 



dr2 



+ + C3if/ + C32(r)6, = 0, (20) 
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where: Cn = P, C22(r) = 1 + (1 + P)K'yiT), C33 = 1 + (1 + P)Kl C^ir) = C2i(r) = 
-f3Ky{r), Ci3 = C31 = and C23(r) = C32(r) = (1 + f3)Ky{T)K,. 

Combining Eqs. (fTH |) - (Pm) it is easy to show that the hnear evolution of the perturbations 
has the following invariant: 

. db* . 



dd Ad* dby d6; d6, 
dr dr ^ dr ^ dr dr 



dr 



tnv, 



(21) 



Here and hereafter asterisk denotes complex conjugated value. Eq. represents conser- 
vation of wave action of the system 0, 0, ■ 

Eqs. ()18 |) -()20 |) describe the linear dynamics of MHD waves in the constant shear flow. 
However, variables d, by, bz are not normal. This circumstance complicates physical treat- 
ment of perturbation dynamics. Therefore, in the next subsection normal variables are 
introduced in the shearless limit. Afterwards Eqs. (fTH|l - ()2()|l are rewritten in the normal 
variables and general features of perturbation dynamics as well as particular characteristics 
of the system evolution for different values of parameters (3, S and are studied in detail. 

A. Shearless limit 



For the sake of clearness of further analysis first we shortly discuss the shearless limit. 
Introducing normal variables \l'j, {i = 1, 2, 3) {2^ : 



(22) 



where 



by 
ybzj 



(23) 



^ cos a sin a cos 7 sin a sin 7 ^ 



V 



Q = sin a — cos a cos 7 — cos a sin 7 . (24) 
sin 7 — cos 7 

a and 7 stand for Euler angles in the ■i/'-space. After quite long but straightforward algebra 
it can be shown, that: 



a = arctan 



^1-/3 



(25) 



7 = arctan ( ) , (26) 



with Kl = Kl + Kl 

In the normal variables Eqs. (fTHjl - ipUjl are decoupled and reduced to 



Vti^i = 0. (27) 



dr2 

These equations describe independent oscillations with so-called fundamental frequencies 



1±4 1 



4/3 



nl = 1, (28) 

that are eigen-frequencies of the set of Eqs. (fTBj) -()20j ) in the shearless limit and can be 
easily identified as fast and slow magnetosonic (FMW and SMW) and Alfven waves (AW), 
respectively. Hereafter parallel with notations of eigen-frequencies Qi and eigenfunctions 
\E'j where i = (1,2,3) we also use the notations associated with corresponding MHD wave 
modes, i.e., hereafter indexes of Qi and \E'j vary over (1, 2, 3) or equivalently {F, S, A). 
Substituting Eqs. (j^ - PHj) into one can get the following expressions for the normal 

' ^ ^{n'i - f3f + f3^Kl ' 

^ _ PKlp+{nl-pK^){Kyby + K,K) 

^3 - = ^i^^^. (31) 
B. Dynamical equations in normal variables 

In shear flows, coefficients in Eqs. (fTHjl - ()20p vary in time. Therefore the equations 
governing the dynamics of SFHs become coupled [in contrast to the shearless equations 
(I27|l ] and take the following form in the normal variables: 

^ + [n] - 0,) VI/, = -T;u^, - A,fc^ , (32) 



where: 



(33) 
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QijQkj i ^ k 

i = A; 



i = A; 



(34) 



(35) 



Here and hereafter overdot denotes r derivative. f2j and Qij are defined by the same expres- 
sions as in the previous subsection [see Eqs. (j21I)-(I2EI) and (j2HI)], but with time dependent 
normahzed wave number Ky{T) = Ky — St. Expressions for the coefficients in Eq. (jH^ are 
given in the Appendix A. 

There are three main differences between Eqs. flS^ and shearless equations (jTTj) caused 
by the time dependence of Qij in the shear flow. These differences correspond to three chan- 
nels of energy exchange processes and are responsible for novel features of linear dynamics 
of the perturbations in the shear flows: 

- time dependence of the eigen-frequencies [VLi = Qi{Ky{T))] causes adiabatic energy ex- 
change between main flow and perturbations; 

- terms on the right hand side of Eq. ()32|) describe the coupling between different wave 
modes; 

- additional terms (Oj) on the left hand side of Eq. (j32|) describe shear induced modiflcation 
of frequencies. Influence of these terms become remarkable at high shear rates {S > 1) and, 
as it will be described in Sec. V, they are responsible for the over-reflection phenomenon of 
MHD wave modes. 

Our further efforts are focused on the analysis of this novelties of the perturbation dy- 
namics. 

C. Adiabatic evolution of SFHs 

There are two necessary conditions that should be satisfled for the validity of WKB 
approximation. Firstly, shear modifled frequencies Clf = Q'^ — 0j must be slowly varying 
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functions: 



(36) 



Second condition implies that coupling terms in Eq. (jH^ have to be negligible. 
Let us first consider the limit 5* ^ 1. In this case condition ()36|) reduces to 



(37) 



This condition indicates that WKB approximation fails in some vicinities of the turning 
points where ^^(r) = 0. But direct evaluation of Eqs. ()28p indicates that none of the turning 
points are located near the real r-axis and therefore the condition ()37|) is satisfied for SFHs 
of all wave modes at any moment of time for arbitrary values of the parameters Ky{T), 
and (3. 

If the frequencies of two oscillating modes became equal at some moment of time, Q 
becomes degenerated. Equivalently, as it can be seen from direct evaluation of Eqs. ()A2|) - 
(IA3|) . coupling coefficients of corresponding oscillations in Eq. ()32j) which otherwise are of 
order S or S"^ become infinity. So, at 5 ^ 1 evolution of SFHs of the wave modes is adiabatic 
except some vicinities of resonant points where f2j(r) = Vljij), or equivalently during some 
time intervals when frequencies of different wave modes are close to each other. Analysis 
of Eq. ()28|) yields that all the resonant points are located on the axis Re[-R'j,(r)] = in 
the complex r-plane. Consequently, these time intervals can appear only in a certain /S.Ky 
vicinity of the point Ky{T) = (exact conditions for effective coupling between different 
modes will be formulated in the next section). 

Now consider the case of moderate and high shear rates when condition S* <^ 1 is not 
satisfied. Again, analysis of turning points where Cli{T) = and resonant points where 
Cli{T) = 0,j{T) is important. But in contrast with the limit 5^1, there is no small param- 
eter in the problem. Therefore, in general WKB approximation may fail during the whole 
evolution. Combining Eqs. ()28|) and ()25 p - (j^Bj) after long but straightforward calculations 
one can conclude, that both conditions for the validity of WKB approximation formulated 
above hold for SFHs of all MHD wave modes for arbitrary (3 and at least for high values 



This circumstance is crucially important for the analysis of the mode coupling at high shear 
rates presented in Sec. V. 



of \Ky{T)\: 



ir,(r)|»5. 



(38) 
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If the conditions for the vahdity of WKB approximation holds, temporal evolution of 
SFHs can be described by standard WKB solutions: 

n± 

where Df are WKB amplitudes of the wave modes with positive and negative phase velocity 
along X-axis respectively. All the physical quantities can be easily found by combining Eqs. 
(j29|) - (j3ip . Combining these equations after long but straightforward algebra one can check, 
that the energies of the wave modes satisfy standard relations of adiabatic evolution: 

E, = n,{T){\Dt\' + \D7\% (40) 

so asymptotically (r —>■ oo) total energy of SFH of SMW and AW become constant whereas 
the total energy of SFH of FMW increases linearly with r. In other words, SFH of FMW 
can effectively extract energy from the basic flow. It follows from Eq. pUj) . that l-D^^P can 
be interpreted as number of wave particles (so called plasmons) in analogy with quantum 
mechanics. 

Using Eq. (jH^ we can reduce Eq. (|7T|) to the following form: 

J2\Dt\'-Y^m' = ^nv. (41) 

i i 

This equation is crucially important for the study of non-adiabatic processes in the consid- 
ered flow: as it was mentioned above, the evolution of SFHs of the wave modes is always 
adiabatic for arbitrary Ky{T) which is outside some AKy vicinity of the point Ky{Tj^) = 0, 
i.e.: 

\KyiT)\>AKy, (42) 

where the value of AKy depends on the specific parameters of the problem. From the point 
of view of temporal evolution this means that non-adiabatic processes can be important only 
during the time interval Ar = AKy/S in the vicinity of r+, where condition ()42|) fails. 

Consequently, the dynamics of SFHs is the following: assume at the initial moment of time 
Ky{0) > AKy. According to Eq. (jl^ the dynamics of SFHs is adiabatic initially. Due to 
the linear drift in the k-space, Ky{T) decreases and when condition fails, the dynamics 
of SFHs become non-adiabatic. Duration of non-adiabatic evolution is Ar. Afterwards, 
when Kylr) < —AKy, the evolution of SFHs becomes adiabatic again. 
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Denote WKB amplitudes of SFHs of the wave modes on the left and right sides of the 
area of non-adiabatic evolution (i.e., for r < r_|_ — Ar/2 and r > r+ + Ar/2) by Df^ and 
D"^^ respectively. In other words, D^^ and Df^^ are WKB amplitudes before and after 
non-adiabatic evolution. Eq. (j41|) provides important relation between WKB amplitudes 
on the left and right sides of the area of non-adiabatic evolution, independent of the be- 
haviour/dynamics of the system in the non-adiabatic area: 

E i^rj' - E \D^,L\' = E \DU\' - E i^^^i'- (43) 

i i i i 

III. GENERAL PROPERTIES OF MODE COUPLING 

As it was mentioned above WKB approximation is valid for arbitrary Ky{T) except some 
vicinity of the point Ky^r^). In the formal analogy with S-matrix of the scattering theory 
2^ and transition matrix from the theory of multi-level quantum systems one can 

connect Dfj^ with Dfj^ by 6 x 6 transition matrix: 



'R 

where and are 1x3 matrices and T^^ are 3x3 matrices. 




(44) 



Due to the 



and — 



act that all components of the matrix C in the set of Eqs. (fTH |) - (Pm) are real 

7;++ = [T,T-]* ^ (45) 



= iTij^T ^ T,,. (46) 
Substituting Eq. (gH) in Eq. dHJ one can get: 

Ei^^.-i'-Ei^^/ = i (47) 

j j 

In general, the components of transition matrix in Eq. (j44p are complex. This means, 
that the interaction of different wave modes changes not only the absolute values of Df, but 
also their phases. We call the absolute value of the transition matrix components \Tij\ and 
\Tij\ the transformation coefficients of corresponding wave modes: 

- \Tij\ represents the transformation coefficient of j to i mode, that has the same sign of the 
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phase velocity along X-axis (i.e., transmitted mode i); 

- \Tij\ represents the transformation coefficient of j to i mode, that has the opposite sign of 
the phase velocity along X-axis (i.e., reflected mode i). 

nn 

It has to be noted that in the similar problems of quantum mechanics [2a,|22|, and 
|Tjjp represent transitions probabilities between different quantum states. 

As it follows from Eqs. ()4Up and ()44|). if initially only one, for instance j mode exists, 
i.e., Df]^ = for 2 7^ j, the energies of transformed waves do not depend on the phases of 
transition matrix elements and are entirely determined by \Tij\ and \Tij\. In the presented 
paper we concentrate attention on transformation coefficients and do not regard the problem 
of the phase multipliers of transition matrix elements. 



Now consider t 
well known l26 



le case 5* <^ 1, for which the condition (j36|) holds at arbitrary r. It is 
27|, that in this case components of T are exponentially small with respect 
to the large parameter 1/S and can be neglected. Consequently, Eq. decomposes and 
reduces to: 

D+ = TD+, (48) 
D« = T*D^. (49) 

From physical point of view this means, that only wave modes with the same sign of the 
phase velocity along X-axis can effectively interact, i.e., the wave reflection is negligible. 
Eq. (jlTj) reduces to the unitary condition for T: 

J2\T^J\' = l. (50) 
j 

Due to the fact that condition is satisfied, the only reason for the failure of WKB 
approximation could be closeness of the wave mode frequencies, i.e., closeness of at least one 
resonant point (point where two of more fundamental frequencies become equal) to the real 



r-axis j26|, |27[ . In the later case coupling of wave modes becomes effective. Analysis of Eq. 
(j28|) shows that all the resonant points are located on the axis Re[ii'j^(r)] = in the complex 
r-plane. That is why the effective transformation of wave modes can take place only in the 
vicinity of the moment of time r+ where Ky{T+) = 0. 

First of all let us discuss some general properties of wave resonant interaction (mathe- 
matical details are discussed in the Appendix B): 
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(i) for effective coupling between different (for instance i and j) wave modes there should 



exist a time interval (so-called resonant interval) where 



26 



(51) 



If this condition is not satisfied transformation coefficients are exponentially small with 
respect to the large parameter 1/S, namely 0, Eh 



Tij ~ 



exp 



Im / {^li — Vtj)dT 



(52) 



Here and hereafter the signs of absolute magnitude for transformation coefficients are omit- 
ted, i.e., under T^j we mean |Tjj|. In Eq. (j52j) . tq is arbitrary point on the real r-axis and 
Tij is the nearest to the real r-axis resonant point where fij(rjj) = f2j(rjj). 

The characteristic equation of the set ()18 |) - (j^m) is real and symmetric with respect to 
the transform Ky{T) —Kyij). Therefore, resonant points always appear as complex 
conjugated pairs: if T^j is a resonant point, so is its complex conjugated one r-*. 

(ii) the following equation holds for resonant interaction of two wave modes, e.g., i and j 
(i.e., when condition ()51|1 is satisfied only for two wave modes): 



T • = T • 



(53) 



This symmetry property, which follows from unitarity property (jSDI), holds for resonant 
interaction of two wave modes only. If at the same interval of time there is effective coupling 
of more then two wave modes, then Eq. fails. 

(iii) if in the neighborhood of the real r-axis only a pair of complex conjugated first order 
resonant points Tij and T*j exists [the resonant point is called of order n if fQ? ~ ^f) ~ 
{r - niT'"^ in the neighborhood of r.,], the transformation coefficients are HQ: 



T- ■ 

7.1 



exp 



Im / {Qi — Vtj)dT 



TO 



\1 + 0{S^'% 



(54) 



There are two remarks about this equation. Firstly, it shows, that only dispersion equation 
is needed to derive the transformation coefficient with accuracy S^^'^ in the case of the 
first order resonant points. In other words, only the solution of the characteristic equation 
of the governing set of Eqs. (f7j)- (fT3j) is needed to derive the transformation coefficient. 
Secondly, Eq. (jMj) is valid also in the case of strong wave interaction. For example, if 
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complex conjugated resonant point of the first order tends to the real r-axis, then T-ij 1. 
According to Eq. (jSUjl . this means that one wave mode is fully transformed into another. 

(iv) for the second or higher order resonant points analytical expression of transformation 
coefficient can be derived only in the case of weak interaction (Tjj ^1, i ^ j). Namely, if 
there exist a pair of complex conjugated second order resonant points and condition (j51|) is 
not satisfied, transformation coefficients are: 

T,, = ^exp Im^ - n,)dr ^ [1 + 0{S'/% (55) 

Let us note once again, that this equation is valid only for Tij -C 1. Derivation of this 
formula is presented in Appendix B. 

In earlier studies, an attention was always paid to the resonant points of the first order 
27l |. As it will be shown later, all the resonant points are of the second order in the 



case of resonant coupling between SFHs of AW and the magneto-sonic wave modes. This is 
not an unique property of the evolution of MHD wave modes. It can be readily shown that 
these type of resonances naturally appear in systems of three or more coupled oscillators. 

(v) from Eq. ()51|) it can be shown (see Appendix B), that the time scale of the resonant 
interaction is as follows: 

An, ~ (56) 

where n is the order of the resonant point. Whereas the time scale of adiabatic evolution 
At ~ 1/^. Therefore if 5 < 1 

ATij < At, (57) 

i.e., resonant interaction of waves is much faster process then energy exchange between 
background fiow and wave modes. Consequently, conservation of wave action (jTrf) reduces 
to energy conservation during the resonant interaction of wave modes. 



IV. SPECIFIC LIMITS OF THE RESONANT TRANSFORMATION OF MHD 
WAVE MODES 

In this section we study specific cases of MHD wave coupling at S* ^ 1: 
(A) Two dimensional (2D) problem when = 6^ = 0. In this case Alfven waves are 
absent, and obtained set of equations describes the coupling between FMW and SMW. 
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(B) f3 <^ 1. In this limit Qf,a S> and only mutual transformation of FMW and AW 
is possible. 

(C) /5 1. In this limit flA,s fip and mutual transformation of SMW and AW can 
be effective. 

(D) /? ~ 1. In this case frequencies of all the MHD waves can be of the same order and 
mutual transformations of all the modes is possible. 

As it was mentioned above resonant transformations of MHD wave modes in shear flows 



are investigated recently ll3|. The content of this section is concerned on detailed 



quantitative analysis of the problem. In particular - derivation of analytical expressions of 
transformation coefficients. 



A. 2D problem 



To derive equations in 2D case one has to assume 
Alfven waves and Eqs. (fTHjl - lpljl reduce to: 



+ Cud + Ci2ir)by 



(j^2 ' ^^^^ ' ^^^v' j-^y 0) 



0. This limit excludes 



(58) 



dr2 



+ C22iT)by + C2l{T)d = 0, 



(59) 



where: Cn = P, C22(r) = 1 + (1 + /?)i^2(r), Ci2(r) = C2i(r) = -{SKyir). Eqs. 
describe coupled evolution of SFHs of FMW and SMW. Frequencies of the modes are given 
by Eq. (j2HI) where now K"^ = 1 + Kyi^r). Normal variables are defined by Eqs. pHjl - ipiljl 
with substitution = = 0. 

Solving the equation Qp = one can easily obtain that there is only a pair of complex 
conjugated resonant points of the first order in the complex r-plane: 



Ky{,TFs) 



(3+1 



^y\' FS) 



—%- 



Noting, that in the neighborhood of resonant points: 



from Eq. one can easily obtain: 

Tps ~ exp 



1 1/2 



4^25 



(60) 



(61) 



(62) 
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It is seen from Eq. fl6U|) . that if /9 ^ 1, the resonant points tends to the real r-axis. Then 
from Eq. (jU^ it follows that Tps 1- According to Eq. (|3n|). it means that one wave mode 
totally transforms into another. 

Dependence of transformation coefficient on (3 is presented in Fig. Q for S* = 0.05 and 
S = 0.02. Dotted line shows transformation coefficient obtained by numerical solution of 
the set of Eqs. (fTHj) -()20 |) . where initial conditions are chosen by WKB solutions (jH^ far 
on the left hand side of resonant time interval {Ky{0) ^ — 1) and solid line is the curve of 
analytical solution 

B. Low l3 regime 

In this case the frequency of SMW is far less then frequencies of FMW and AW {Qs ^ 
Qf,Qa)- Therefore, the coupling of SMW with other MHD modes is exponentially small 
with respect to the parameter 1/S and can be neglected. Consequently, in the set of Eqs. 
(fT^ - pnjl equation for d decouples and equations for by and describes coupled evolution 
of SFHs of AW and FMW: 

^ + (1 + K'yir)) by = -Ky{T)K^b., (63) 

+ (1 + K',) b, = -Ky{T)K,by, (64) 
Normalized frequencies of the coupled wave modes are: 

nUr) = l + K',+K'^{r), ^1 = 1. (65) 

From this equation it follows that there are two complex conjugated second order resonant 
points: 

Ky{TFA) = iK„ Kyir*^) = -iK,. (66) 
Necessary condition for effective coupling expressed by Eq. 1)511) now takes the form: 

\Kl\< S. (67) 

Thus, the critical parameter is 5 = Kz/S^^^. 
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Consider the case 5 3> 1. For the calculation of the transformation coefficients one can 
use general analysis presented in Sec. Ill and Appendix B. Specifically, Eq. yields: 

T^A=2exp^ ^ J' (68) 

where: 

(t>FA{K,) = (1 + Kl) arctan(ir,) - K,. (69) 
If in addition Kz <^1, then Eq. reduces to: 

T^A-|exp(-0. (70) 

Analytical expression for the transformation coefficients can be derived also in the oppo- 
site limit (5^1, {Kz ^ S^^'^). Taking into account the relation fi| /?, one can readily 
obtain from Eqs. and (j^Tj) that hy and hz coincide with the eigenfunctions of SFHs of 
FMW and AW respectively accurate to the terms of order K^. Consequently, terms on the 
right hand sides of Eqs. represent the coupling terms with accuracy K^. Since 

Kz <^ 5*^/^, coupling is weak and feedback in the set of Eqs. (j63|) - (j6H) strongly depends on 
the amplitudes of the wave modes. Then it follows that if initially only AW exists, feedback 
of FMW on AW is negligible and the set of Eqs. reduces to: 



dr2 



+ [1 + irj(r)] = -Ky{T)Kz^A. (71) 



^ + = 0. (72) 



dr2 

Using the solution of Eq. ()72j) and well known expressions for the solution of linear 
inhomogeneous second order differential equation, in the considered limit (5 ^ 1), we obtain: 

Tfa ~ j X sin {- - J dx. (73) 

Note that: 







X sin ( — 7xJ dx = tt— y4z(— 7) (74) 



and using the expansion of Airy function Ai^'y) into power series [29| we finally obtain: 

2^/^- J. r(|) 



-FA 



3V3r(i)M^ 2V43i/3r(|)'^ '■ (^^^ 
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Results of numerical solution of the initial set of equations ff7|)- (fT^ (solid line) as well as 
analytical expressions (fTOj) (dash-dotted line) and (f73j) (dashed line) are presented in Fig. 
121 It shows that the transformation coefficient reaches its maximal value {Tp^)max = 1/2 at 
6'^'^ that can be found numerically or alternatively by finding the maximum of the analytical 
expression presented by Eq. (f7H|) : 

1/4 

(76) 



'27/43l/3rgy 



5r(|) 

Eq. (f7F)|l is in perfect accordance with numerically calculated 5"^^ (see Fig. ^ despite the 
failure of Eq. ()75|) at 5 ~ 1. This fact can be explained as follows: the only reason of failure 
of Eq. (f73)) is the neglect of the feedback in Eq. ((721). The feedback changes the value of 
the transformation coefficient but does not affect on the value of 6'^^. 

{TFA^)max = 1/2 mcaus that only half of the energy of FMW can be transformed into 
AW and vice versa even in the optimal regime. It has to be noted, that Landau-Zener 
theory jl?! . Q] provides the same maximum value for the transition probability in two-level 
quantum mechanical systems. 



C. High (3 regime 

In this case fi^, VLa ^^f- Consequently, the coupling of AW and SMW with FMW are 
exponentially small with respect to the parameter 1/S and can be neglected. 

Analysis of Eq. ()28|) provides that for the coupling of AW and SMW there exist two 
complex conjugated second order resonant points (solutions of the equation VLs = ^a), that 
are also given by Eq. The condition of effective coupling (j5ip takes the form: 

\K P 

<f3S. (77) 



If this condition fails, then Eq. ()55|) yields exponentially small transformation coefficient: 

TsA = - exp — , (78) 



2I3S 

where 

arcsinh [K^ 



sa{K.) = K^- r-—^ - (79) 
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If additionally Kz ^ I, then Eq.(f7Hj) reduces to: 



7r 



T5A^-exp(-i^). (80) 

Consider the transformation process when the condition of effective coupling Eq. (ffTjl 
is satisfied. Similar to the low (3 regime, in this case it is also possible to derive the set 
of two second order coupled equations that describe the linear coupling of AW and SMW. 
Expressing the density perturbation p from the condition = (thus eliminating the 
terms that describe interaction of AW and SMW with FMW) and combining Eqs. (|T8|) - (j2n| 
one can obtain: 

^ + (1 + OK'^ir)) by = -eKy{T)KA, (81) 
^ + (1 + ^i^f) K = -eKy{T)K,hy, (82) 



where: 



As it follows from Eqs. ()81|) - (jH^ Tas = at = 0, i.e., there is no interaction between 
AW and SMW. 

In contrast with low /3 limit, there is no unique parameter in high f3 limit that totally 
describes the transformation process. In this limit there are two such parameters and 
f3S. 

First consider the case pS <t:l. Then the condition ^ reduces to 6i = \K,\/{pSy/^ < 
1. Note, that the sign of 6 does not affect the transformation coefficient and K'^{r) ^ 1 in 
the resonant area. Thus we conclude that the properties of the wave transformation is the 
same as in the case of transformation of AW and FMW. Namely, if 5i <^ 1, then the leading 
terms of asymptotic expressions of transformation coefficient is given by Eq. (|75p with S 
replaced by 6i: 

At (3S <^ 1, the transformation coefficient reaches its maximum {T'^s)max = 1/2 at 6f 
that coincides with 5"' defined by Eq. (fTB)) . 

Dependence of the transformation coefficient on Kz obtained by numerical solution of the 
initial set of equations (fTH j) -(p? j) for pS = 0.025, pS = 0.5 and pS = 1 are presented in Fig. 

m 
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Consider the case when /3S is not small. The numerical study shows that in this case the 
properties of transformation process is totally different (see Fig. 3): 

- transformation coefficient Tas does not depend on jSS at i^r;^ -C 1 

Tas ~ 2.05ir„ (85) 

- (TAs)max = 1, i-G., total transformation of one wave mode into another is possible. 

D. /? ~ 1 regime 

In the case of /3 ~ 1, the frequencies of all MHD wave modes are of the same order and 
no simplification of the set of Eqs. (fTH jl -fpU j) is possible. Analysis of Eq. shows that 
there exist the pair of complex conjugated first order resonant points: 

Ky{n,2) = ±^^i^|+(|^y (86) 
and the pair of complex conjugated second order resonant points: 

Ky{rs,4) = ±tK,. (87) 

No analytical expressions can be obtained for transformation coefficients if more than 
two wave modes are effectively coupled. 

Numerical study of the set of Eqs. (fTHjl - lj^Uj) is performed as follows: WKB solutions 
()39|) are used to obtain initial values of d, by, and their first derivatives for different wave 
modes separately far on the left side of the resonant interval {Ky{0) ^ 1). After passing 
through the resonant interval (i.e., for any r, for which Kyij) ^ — 1), WKB solutions were 
used again to determine the intensities of the transformed wave modes. 

Figure m shows results obtained by numerical solution of the set of Eqs. flTH jl - ipUj) . for 
initial SMW. Namely, transformation coefficients (Tps,Tas,Tss) vs are presented for 
different values of /3 and S. 

According to our numerical study, qualitative character of wave transformation process 
is similar to the cases described in the previous sections. However, there are some differ- 
ences. The most interesting is the failure of symmetry property (jSSI)- The presence of third 
effectively interacting wave mode leads to the asymmetry of two wave mode interactions, 
e.g., intensity of AW generated by SMW differs from the intensity of the inverse process. 
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V. HIGH SHEAR REGIME: WAVE OVER-REFLECTION 



Let us consider the dynamics of SFHs of MHD wave modes in the flow with moderate and 
high shear rates (5* > 1). In this case the dynamics is strongly non-adiabatic [Gj terms can 
not be neglected in the left hand side of Eqs. ((221)] except time intervals when |ii"y(r)| ^ S. 
The existence of these adiabatic intervals permits to study wave interaction based on the 
asymptotic analysis presented in Sec III. 

The main novelty that appears in the flow at high shear rates is that the dynamics 
involves wave reflection/over- reflection phenomena. Specifically, Tij can not be neglected 
and becomes important in Eq. (j77|) . Physically it means that initial "^f mode can be 
effectively transformed into "^J modes (wave modes with phase velocity directed opposite to 
phase velocity of initial wave mode with respect to X-axis). Thus, at high shear rates, the 
wave dynamics represents an interplay of transformation and (over)reflection phenomena 
that are described by Tij and Tij. 

Dependence of T^j and Tij on obtained by numerical solution of the initial set of 
equations (fTS|) - (f2m) for different values of parameters (3 and 5* are presented in Figs. EHHl 
Figs. EEshow the cases when initially only SMW with positive phase velocity along X-axis 
exists, i.e., only D^^ ^ 0. Whereas in Fig. IHlinitially only AW with positive phase velocity 
along X-axis exists, i.e., only D\j^ ^ 0. 

In Fig. = 2 and S = 1 (dashed lines) and S" = 4 (solid lines). Both, wave 
transformation and reflection phenomena appear and enhance with increase of S. Coupling 
of SMW with FMW is dominant at i^^ -C 1, while coupling of SMW with AW is dominant 
at Kz > 1. Figure El shows the same plots for high (3 and moderate S. Particularly, dashed 
lines correspond to /3 = 50 and 5* = 1, whereas solid lines correspond to /3 = 50 and 
S* = 2. For this range of parameters ^ S and according to Eq. ((U^ Tps^Tps ~ 0. 
Hence, SMW and FMW are not coupled. Fig. IHl shows one more new feature of the wave 
dynamics. Namely, Tss > 1 at 5* = 2 and Kz <^ 1. Physically it means that the amplitude 
of the reflected SMW is greater than amplitude of incident SMW. This is the over-reflection 
phenomenon first discovered by Miles for acoustic waves [s^. This phenomenon becomes 
dominant at S* ^ 1. 

Figures [7| and |H1 show the case when (3 = 100 and S = 8, 12, 16, 20. The dashed lines 
correspond to S* = 8. For this range of the parameters all the wave modes are coupled. The 
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wave over-reflection phenomenon is more profound. Particularly, in Fig. n\Tss,Tss ~ 7 
when Kz <^ 1. In other words, the energy density of the reflected (as well as transmitted) 
SMW is about 50 times greater that the energy density of the incident SMW. When Kz ~ 1, 
then Tas,Tas ~ 5, i.e., the reflected and transmitted AW are substantially greater then 
incident SMW. Figure |H1 represents the transformation coefficients if only AW exists initially. 
In this case significant growth of energy density of perturbations takes place at Kz > 1. 

For more detailed analysis of the over-reflection phenomena consider the dynamics of 
SMW at Kz = 0, /3 > 1 and 5 < p^^^. Due to the condition Kz = there is no coupling 
between SMW and AW. Condition S <^ /S^^^ provides that the coupling between SMW 
and FMW is also negligible. Combining Eqs. (jH^ . and ()A1|) - ()A3|) one can obtain following 
second order equation for the evolution of SMW: 



dr2 

where: 



f{S.T)<Ss = 0. (88) 
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Asymptotic analysis of this kind of equations is well known in quantum mechanics 
and the theory of differential equations 13- From mathematical point of view the 



reflection is caused due to the cavity of f{S, r) that appears in the vicinity of the point 
where Ky^r^) = 0. 

In accordance with quantum mechanics we define reflection and transmission coefficients 

-2 n n n 

as -R = Tgg and T = respectively. It is well known |26l . l27l |3l| that wave action 
conservation provides that: 

1 = Tis - ^T-R. (90) 

So we can conclude that the total energy of SMW always increase during the reflection/over- 
reflection process: 

^ = T+R. (91) 

The dependence of reflection coefficient Ron S at (3 = 400 obtained by numerical solution 
of Eqs. ()18|) - (j^n|) is presented in Fig. El Initial conditions for numerical solution are chosen 
by WKB solutions far on the left side of resonant area {Ky{0) ^ 1,5'). For high 
shear rates {S ~ P^^'^) SMW is partially transformed to FMW, that is why the growth rate 
of the reflection coefficient decreases (dashed part of the graph in Fig. According to 
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the numerical study of Eq. (|HHjl. the amphtude of reflected SFH exceeds the amplitude of 
incident SFH if: 



This condition indicates that the phenomenon of over-reflection can take place in plasma 
with /? 3> 1 even for small values of the shear parameter A. The amplification of the 
energy density of perturbations is always finite, but it can become arbitrarily large under 
due increase of the shear rate. 

VI. SUMMARY AND DISCUSSIONS 

The linear dynamics of MHD wave modes in constant shear flows is studied qualitatively 
and quantitatively in the framework of non-modal approach. Usage of asymptotic analysis 
familiar in quantum mechanics allows to study the physics of the (over) reflection and mutual 
transformation phenomena of MHD wave modes. Quantitative asymptotic and numerical 
analysis is performed for wide range of the system parameters: relation of spanwise and 
streamwise wave numbers {Kz = kz/k^), plasma beta {(3 = c^/V^) and the shear rate 
{S = A/{kxVA))- The transformation takes place at small as well as at high shear rates 
and involves all MHD wave modes. The over-reflection becomes apparent only for slow 
magnetosonic and Alfven wave modes at high shear rates. 

At 5 ^ 1, the transformation has resonant nature and different wave modes are coupled 
at different (3. 

- In high f3 regime Alfven and slow magnetosonic waves are coupled and the transformation 
is effective at — 1. At (33 <^ 1, for the maximum transformation coefficient we have 
(T^s)max = 1/2. If PS is not small, {TAs)max = 1, i-e., the total transformation of one wave 
mode into another is possible. 

- In low j3 regime Alfven and fast magnetosonic waves are coupled. The transformation is 
effective also at — 1. (TpA^)max = 1/2, i.e., even in the optimal regime only half of the 
energy of SFH of FMW can be transformed into AW and vice versa. 

- At /5 ~ 1 all of the MHD wave modes are coupled: at -C 1 mainly are coupled slow 
and fast magnetosonic waves. Whereas at Kz — 1 the coupling of these waves with Alfven 
wave is dominant. 

At moderate and high shear rates (S), the wave mutual transformation is accompanied by 




(92) 
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wave (over) reflection phenomenon. The amphfication of the energy density of perturbations 
is always finite. However, the ratio of final to initial energies of waves can become arbitrarily 
large with increase of the shear rate. 

Described phenomena could have substantial applications in the theory of MHD turbu- 
lence that is intensively developing recently (see, e.g., Ref. js^ and references therein). This 
concerns both the fundamental problem of MHD turbulence and wave modes participation 
in the turbulent processes. 

Fundamental problem is the elaboration of the concept of self-sustaining MHD turbulence. 
None of the previous studies were concerned on energy source of the turbulence. The usual 



procedure is the following [32]: some energy source of perturbations is assumed a priori 
and the evolution of energy spectrum is studied. On the other hand, in the framework 
of canonical MHD theory, there is no linear instability in the flows (they are spectrally 
stable). Therefore, existence of an energetic source for the permanent forcing of turbulence 
is problematic in the canonical/spectral theory. Nevertheless, one can apply the concept of 
bypass turbulence elaborated by hydrodynamic community for spectrally stable shear flows. 
The role of the energy source in the MHD flows could play waves non-exponential growth 
(over-reflection) . 

Due to the importance of over-reflection phenomenon as a possible energy source for 
MHD turbulence in spectrally stable flows, detailed analysis of evolution of SFHs of MHD 
wave modes in flows with high shear rate will be treated more extensively elsewhere. 

As for ingredients of the turbulence: the study of homogeneous MHD turbulence provides 
that in the case of /5 ^ 1, the turbulence is mainly Alfvenic (i.e., mainly consists of AW) 
and SMW plays a role of passive admixture. Studied in this paper linear coupling of MHD 
wave modes can significantly change this scenario even in the presence of small velocity shear 
{S ^ 1). Thus for the wide range of the system parameters all of the wave modes (at list 
two of them) should be the ingredients of the real MHD turbulence - the turbulence should 
be of a "mixed" type. 
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APPENDIX A: COUPLING COEFFICIENTS 



Substituting Eq. 
obtain: 



into Eqs. after straightforward calculations one can 



^ + 7^ sin^ a ^ 







dp' + 7^ cos^ a 



\ 



r 



/ 



A 



v 



a 7 sin a 

— d —7 cos a 
-7 sin a 7 cos a 



and for non-zero components of T: 

T12 = d^ + 7 sin a cos a, 

Ti3 = 72sina, 

T21 = — d + 7^ sin a cos a, 

^^23 = -7 cos a, 

T31 = — 2d7cosa — 7sina, 

T32 = — 2d7sina + 7cosa, 
where a(r) and 7(r) are defined by Eqs. and with Kyij-) = Ky — St. 



(Al) 



(A2) 



(A3) 



APPENDIX B: SOME MATHEMATICAL ASPECTS OF LINEAR TRANSFOR- 
MATIONS 



To avoid the complication of mathematical formalism first consider the case of two coupled 
oscillators: 



dV 



dr2 

dr2 

with G12 = G21 and eigen-frequencies: 

1 
2 



r + Gii^i = Gi2^/'2, 



+ 022^^2 = G'2lV'l 



^'1,2 



Gn + G22± V (Gil -6-22)2 + 4G' 



12 



(Bl) 



(B2) 



(B3) 
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that are slowly varying functions of r: 

ni<^nl i = i,2. (B4) 

T is normalized so that Qii^Sr) ~ 1 and S <^ 1. In addition, let us assume that all eigen- 
frequencies are real during the evolution. It can be easily seen from Eqs. ()B1|) - ()B2|) that 
this condition implies: 

GnG22 > Gl,. (B5) 

From this condition it follows that neither turning points (where Qi = 0) nor resonant points 
{Ql = Q2) can be located on the real r axis. Transformation matrix to the normal variables, 
similar to Eq. ()24|1 . is: 

I L^wo /\ 0111 /\ 1 

Q = I I , (B6) 

with: 

A = arctan ( ^^^^ ) . (B7) 
In the normal variables the set of equations ()B1|) - ()B2|) reduces to: 

^'i + (nj - A^) ^1 = -2A^2 - A^2, (B8) 
^2+ (^^2-^^)^2 = 2A^i + A^i, (B9) 

where: 

• _ {nl - Gu)dGu/dr - G^^dj^l - Gu)/dT 

According to Eq. ()B3|) condition = On leads to G12 = 0. Consequently, as it was 
mentioned above, only singular points of coupling coefficient are resonant points, where 

Let and be linear independent solutions of corresponding homogeneous equations 
(Eqs. ()B8|) - ()B9|) with right hand sides equal to zero), normalized such that asymptotically 
they converge to WKB solutions: 

^^^2 ~ ^^exp (±t [ f]i,2rfr ) . (Bll) 



a/^1,2 \ J 

The general solution of the set of Eqs. ()B8|) - ()B9|) is {i = 1,2 



^, = Dl^t + Di^;- (B12) 
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For simplicity let us assume, that far on the left side of the resonant area only one, for 
example mode exists, i.e., D2 = 1 and all other WKB amplitudes are zero. Consider- 
ing terms on the right hand side of Eqs. ()B8|) - ()B9|) as external source, using well known 
expressions for the solution of linear inhomogeneous second order differential equation 
and integrating by parts one can estimate the amplitude of generated wave as: 



dr dr 
Necessary condition for effective coupling implies: 

Df ~ 1. 



dr. 



(B13) 



(B14) 



Presenting and in the form of formal series jSlJ, substituting them into Eq. ()B13|) 
and taking into account that obtained series of integrals converges rapidly jsS, 34|, condition 
flB14|) takes the form: 



A exp 



(B15) 



If S* ^ 1, then phase multiplier in ()B15|) is rapidly oscillating. Condition ()B14|) can be 
satisfied only if there exists an interval where time scale of the changing A and the phase 
multiplier in the integrand of Eq. ()B15|) are at least the same. This lead to the condition: 



(B16) 



at some interval along the path of integration. On the other hand one can note that only 
singular points of A in the complex r-plane are resonant points (where Qi = ^2)- If so. 
Landau's rule 3| provides that if S* ^ 1, then leading term of asymptotics of transformation 
coefficient is: 



12 



exp 



Im 



(B17) 



{ill - VL2)dt 

where tq is arbitrary point on the real r-axis and T12 is the nearest to the real r-axis resonant 
point. Eq. ()B17|) represents an alternative necessary condition for effective coupling of wave 
modes. Namely, resonant point T12 must be close to the real r-axis in such a way that 
the exponent in right hand side of Eq. ()B17|) has to be of the order of unity. Surely this 
condition is equivalent to Eq. ()B16|) . On the other hand, it shows that if the condition of 
effective coupling is not satisfied, then the transformation coefficient is exponentially small 
with respect to the large parameter 1/S. 
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According to Eq. ()B10|1 . in the neighborhood of nth order resonant point: 

Eq. ()B16|) yields following estimation for the effective coupling time scale: 

Ari2 ~ 5-^+1/". (B19) 

1. Resonant points of the first order 

Consider the case when only a pair of complex conjugated first order resonant points 
exists close to the real r-axis. This means that in the neighborhood of T12: 

\nl - nil {t - TuY^^ (B20) 

In this case, the leading term of asymptotic of transformation coefficient can be found exactly 
by the method first used by Landau 



Consider Eqs. (IB8j) - ()B9j) along the path ( that turns over tu in the complex r-plane 
(Fig fTTH) . Considering Eqs. ()B8jl - ()B9|l in the neighborhood of T12, where A ~ (r — ri2)^^/^, 
one can obtain that Eqs. ()B8|) and ()B9|) exchange after turning over resonant point. This 
means that \I^2 became and vise versa. After this kind of analysis it is usually concluded 
(see e.g., Ref. ^^^^ transformation coefficient is: 

T12 = exp Im^ '\n^ - n2)dT ^ [1 + 0(5^/2)]. (B21) 

From our point of view this kind of analysis needs to be commented. When one considers 
the set of Eqs. flB8|) - ()B9|) along the path (, except mentioned above exchange of and \E'2, 
there are also terms 



c 



dr dr 

In the case of the first order resonant points it can be shown that according Watson's lemma 
contribution of these terms are of order S*"*^/^ exp(— | J^^'^(^li — ^l2)dr\) and ()B21|) is valid. 
However, in the case of the second or higher order resonant points, the contribution of terms 
like ()B22|) is at least of the same order as ()B21|) . That is why this asymptotic method works 
only in the case of first order resonant points. 
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2. Resonant points of the second order 



Unfortunately the transformation coefficient can not be obtained analytically in the case 
of the second or higher order resonant points. Exact asymptotic expressions can be obtained 
only in the case of weak interaction. Let us consider the case when T12 and complex 
conjugated second order resonant points. Assume again, that far away on the left side of 
the resonant area only one, say, \t'^ mode exists. In the case of a weak interaction feedback 
of to Eqs. ()B8|) - ()B9|) can be neglected, so these equations decouple and there 

remains only inhomogeneous second order differential equation. We skip some details of 
calculations noting only main steps. In the case of neglible feedback, for the transformation 
coefficient one can obtain: 



12 



2.1 ^ \ ' dr ^ dr 



(B23) 



In the neighborhood of the resonant point one can use the formal asymptotic series for "^^^2^ 
considering Qi^2{'t) = ^1,2 (T12) as constants. Consequently, the last equation takes the form: 



T12 = ^" y XnA"exp \^ J {Qi - Vt2)dt ) dr 
where Xn are regular functions in the neighborhood of T12 and: 

Xi 



(B24) 



(B25) 



To evaluate the integrals in ()B24|1 . one can use Van der Wearden's method [? ]. It is not 
difficult to show that each term in Eq. ()B24jl are of order 5*^/^ with respect to the previous 
term. Consequently, in the case of the second order resonant point, leading term of the 
asymptotics is defined by the first term, for which we have: 



12 



TT 



im f - n2)dT J . 

•J TO J 



(B26) 
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Figure captions 

FIG.^ Transformation coefficient Tps vs [5 for 5 = 0.05 and S = 0.02 obtained by numerical 
(dotted line) and analytical solutions (solid line), see Eq. (|62j) . 

FIG. [21 Transformation coefficient Taf vs 6. Dash-dotted line and dashed line represent ana- 
lytical expressions (jTOf) and (|75() . respectively. Solid line is obtained by numerical solution of 
Eqs. 

FIG. Transformation coefficient Tas vs for S = 0.005 and = 5 {f3S = 0.025), (3 = 
100 {I3S = 0.5) and /? = 200 {(33 = 1). 

FIG. |1] Transformation coefficients Tps, and T55 vs Kz for /? = 0.8 (left column), [3 = 
1.0 (center column) and (3 = 1.2 (right column) for different normalized shear parameters: 
S = 0.005 (dash-dotted hues), S = 0.02 (dashed hues) and S = 0.08 (sohd hues). Note, that 
according to the conservation of wave action Tpg + T'^g + ~ ^ (^^^ Eg. HT|) . 
FIG. 13 Transformation coefficients {Tis and Tis respectively) vs at /3 = 2 at different shear 
rates: S = 1 (dashed lines) and 5 = 4 (solid lines). 

FIG. El Transformation coefficients Tas,Tas, Tss and Tss vs for /? = 50 at different shear 
rates: S = 1.0 (dotted lines) and S = 2.0 (solid lines). Tps,Tfs ^ 0- 

FIG. d Transformation coefficients {Tis and Tis respectively) vs at /? = 100 at different 
shear rates: 5 = 8, 12, 16, 20. Dashed lines correspond to 5 = 8. 

FIG. IHl Transformation coefficients {TiA and TiA respectively) vs Kz at /3 = 100 and different 
shear rates: 5 = 8, 12, 16, 20. Dashed lines correspond to 5 = 8. 

FIG. ini Reflection coefficient of SMW R{S) vs normalized shear parameter 5. Here P = 400. 
For high shear parameters (5 ~ P^^^) SMW is partially conversed to FMW, that is why that 
slope of the refraction coefficient graph decreases (dashed part of the graph). 
FIG. rrUl Resonant points and path ^ in the complex r plane. 
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